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We study motion of a phase transition (PT) front at a constant temperature between stable 
and metastable states in fluids with the Van der Waals equation of state. We focus on a case of 
relatively large metastability and low viscosity, when no steadily moving PT front exists. Simulating 
the one-dimensional hydrodynamic equations, we find that the PT front generates acoustic shocks 
in forward and backward directions. Through this mechanism, the nonsteady PT front drops its 
velocity and eventually stops. The shock wave may shuttle between the PT front and the system's 
edge, rarefaction waves appearing in the shuttle process. If the viscosity is below a certain threshold, 
an instability sets in, driving the system into a turbulent state. 



Motion of phase transition (PT) fronts is the most important phenomenon in macroscopic physical kinetics jj],^]]. 
In most cases, it is described in the quasi-equilibrium approximation assuming a critical value of the temperature at 
the phase interface 0. A distinct class of phenomena are fast propagating phase transformations, such as detonation 
[3| and waves supported by mechanisms as different as relaxation in glassy solids Q, Belousov-Zhabotinsky reaction 
[5|, bubble boiling 0, and the Marangoni flow coupled to a reaction on a fluid surface or chemical transformations 
coupled to mechanical stress |7) . Coupling of the "chemistry" to the flow or elasticity is a necessary ingredient of the 
fast propagation of a PT wave. 

A related problem is the propagation of condensation and evaporation waves in fluids. In this case, a consistent 
description is possible, based on the combination of the corresponding Navier-Stokes (NS) and heat-transport equations 
with the equation of state, e.g., the Van der Waals (VdW) one |],|l(J. Close to the critical point (CP), the VdW 
equation of a spatially inhomogeneous state takes the following form in terms of the normalized pressure p = (P — 
-Pcp)/-Pcp, temperature r = (T C p - T)/T C p, and density 9 = (p- pep VV^Pcp [§: 

p + At + 3r 3/2 (29 - 9 3 /2) + Gy/¥A6 = , (1) 

where A is the Laplacian, and G is an effective surface tension of the interface separating the vapor (9 < 0) and 
condensed (9 > 0) phases existing below CP (at r > 0). 

The hydrodynamic part of the problem is described by a system of the NS and continuity equations. In the simplest 
case when the PT front is one-dimensional (ID), in terms of the rescaled time t = t-J P c /p c G t (t is the time proper), 
coordinate z = y/r/G x, and flow velocity w = r _1 p c / P c v, these equations are Q 

vh + V^ww z = 69 z - [9/2)9 2 9 z + 9 ZZZ + gw zz , (2) 
9 t + + =0, (3) 
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where g = [(4/3)77 + C] {GPcpPcp) , V and ( being shear and bulk viscosities. 

To produce a closed system, one should add the heat transfer equation to these equations. Analysis performed in 
Ref. |9) demonstrates that a flat PT wave in a bulk medium is driven by the heat transfer giving rise to a slowly 
advancing front. The situation is drastically different in ID case, corresponding, e.g., to the fluid inside a capillar, 
which allows to impose isothermal conditions, r = const. However, it is not necessary to conduct the experiment in 
a capillar; instead, a PT wave may spontaneously propagate in a thin film sheathing an isothermal wire in a bulk 
medium. In fact, a fast PT front of this kind was already observed in experiment ||. 

In the ID situation, a solution for a PT front moving at a velocity —V can be sought for as w = w(£), 9 = 9(£), 
with £ = z + Vt. Eliminating w and assuming 9' (±00) = 0, one arrives at an equation 

9" - gV [1 + v 7 ^ (0 2 - 29)] 9' = (3/2) (9 3 - 9 3 ) 

-[6 + V 2 (l + 2^9 2 )] (0-0 ) , (4) 
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where 9q = 6(—oo) < is the density of the supercooled vapor invaded by the condensation front, and 62 = 6 (+00) > 
is the density of the stable liquid phase behind it; recall that the metastable states have 2/v / 3 < \9\ < 2, while the 
regions \9\ > 2 and \6\ < 2/\/3 are absolutely stable and unstable, respectively. 

An exact solution to Eq. (0) satisfying the boundary conditions 9(±oo) = 6> 2 ,o was found in M, neglecting the 
terms ~ yfr. As is known 0, an equation of this type selects the velocity as an eigenvalue, V(9q). A final result is 
that there is a pair of stable and unstable branches of V(9q), which meet and disappear at a critical density 6* — 
-V^ 8 / 3 ) ( 9 ~ 2 5 2 ) / (6 - 9 2 )- At this point V(6 ) attains its maximum, Vw = I8.9 2 [(9 - 2g 2 ) (6 - g 2 )] ~\ The den- 
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sity of the stable phase behind the front is also found as a part of the solution, #2 = —9^/2+ [(2/3)V 2 + 4 — (3/4)0q] 
In particular, in the case g = (inviscid fluid) 9* = —2, which is exactly the border between the stable and metastable 
states (binodal), hence a steadily moving PT front cannot exist in inviscid fluids. 

Eq. ([|) also gives rise to traveling-wave solutions of a different type, which are acoustic shock waves (ASW) without 
PT. Unlike the PT wave for which the density 6*2 behind the front is determined by 9q, in the case of ASW the densities 
9i and 9f before and after the front may be arbitrary, the velocity being determined by the boundary conditions. 
Neglecting the \fr terms, it is 

Vl sw = (3/2)(9 2 + 9 2 + 9 l 9 f -4) . (5) 

In the limit 8f — $i — ► 0, this yields the sound velocity at a given density, c 2 = (9/2)8 2 — 6. The PT front may 
move faster than sound velocity in the metastable phase invaded by the front. In particular, the velocity V max is 
transonic, provided that g > y/3. Thus, the PT wave is related to detonation, whose characteristic feature is transonic 
propagation [§||. 

The critical density 6** belongs to the metastable region if g < 2. In this case, there is a region of the metastable 
states, 2/\/3 < \9q\ < \9*\, which cannot be converted into an absolutely stable state by a steadily moving PT wave 
||. The objective of the present work is to investigate the propagating PT front in this region by means of direct 
simulations of Eqs. (|^) and (^) (analytical investigation of weakly nonsteady fronts is difficult, as the corresponding 
perturbative expansion contains terms secular in t ||). 

Before that, we notice that, because of the term ~ 99' , Eq. (Q) with yfr terms kept in it does not belong to the 
standard type for which the velocity can be found exactly Nevertheless, an ansatz similar to that solving the 
standard equation yields an exact result for Eq. (0) as well: V = —(3y/3/2g)6i [(1 + \/t9-i) (I + y/r{W\ + 62))] 
9\ 7^ #0,2 being another root of r.h.s. of Eq. ([!]). As yfr <C 1, we look for lowest-order corrections to the final results. 
They contain a new effect, breaking the symmetry between the condensation (+) and evaporation (— ) waves: the 
A/r-corrected threshold value of the viscosity g, below which there is a region of the nonexistence of steady PT waves, 
is 2 ± (5/2)y/r/3. 

In systematic numerical simulations of Eqs. (|2|) and (|3|), the -y/r-corrections were omitted, and w was eliminated, 
differentiating Eq. (|^) in z and substituting w z by — dt as per Eq. (J3J) (for some typical cases, it was checked that 
the yfr terms produce a small insignificant change of the numerical solution). The resulting equation of the second 
order in t and fourth order in z was numerically integrated with the initial conditions 9t(t = 0) = 0, 9{t = 0, z) being 
a smoothed step between the values 9o < and 92 > 0. For the boundary conditions (b.c), it was always taken that 
9 = 9q at the left edge and 9 Z = at both edges. Note that the natural b.c. for the flow, w(— 00) = 0, is compatible 
with these, as is seen from Eq. (|^). The b.c. for the flow at the right edge is not fixed, because, as a matter of fact, 
we are dealing with an open system, with a possible influx of fluid from +00. 

In the case when the steady PT front exists, the fourth b.c. was naturally taken as 9(+oo) = 9 a . Simulations 
have demonstrated that the steady front is always stable when it exists. When it does not exist, the density to be 
established past the PT wave is unknown, therefore initial 92 was taken arbitrarily from the absolute stability region 
\9\ > 2. In this case, the missing b.c. was taken as 9 ZZ = (i.e., the density profile must be flat) at the right edge. 
The results displayed below clearly show that no artifacts are generated by b.c. 

A representative picture of the nonsteady propagation of the condensation front, when the steady one does not 
exist, is shown in Fig. 1. Essential features revealed by many simulations are well seen in it: (i) the nonsteady PT 
wave immediately begins to generate two ASWs in forward and backward directions; (ii) the forward ASW {precursor) 
travels essentially faster than the PT front, so that a rapidly expanding trough with an increased value \9\ > \6q\ is 
formed ahead of the front; (iii) the fast backward-propagating ASW hits the right edge, bounces from it, and then 
hits the PT front. Velocities of all the observed ASWs were checked to comply with Eq. (||). 

Thus, the first finding is that the nonsteady PT wave does not self- accelerate, which could be natural to expect M; 
instead, it generates a strong acoustic precursor. Fig. 2a demonstrates that in all the cases with g < V3 the newly 
established value 9 of the density ahead of the front originally turns out to be fairly close to 0* , corresponding to the 
minimum of \6q\ at which the steady PT wave exists. Then, quasi-steady propagation of the PT wave with the density 
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6 w 9* ahead of it becomes possible. Note that simulations of 2D hydrodynamic equations with the VdW equation 
of state, in the problem of transverse instability of shock waves, have also revealed strong emission of acoustic waves 




The PT front decelerates and eventually comes to halt, while the precursor keeps propagating at a nearly constant 
velocity until it hits the left edge of the integration domain (Fig. 2b). Moreover, simulations show that the PT 
front begins then to slowly move in the reverse direction. We do not consider the reverse motion in detail, as it is 
a manifestation of finiteness of the system, while in this work we focus on properties of a scmiinfinite (to the left) 
system, although, of course, finiteness effects may be important for experiment. As concerns the value of \9\ ahead 
of the PT front, at a later stage of the evolution it slowly increases from the originally established value » |#*| to 
\9\ = 2 at the halt stage. This can be easily understood, as \9\ = 2 is the binodal (a minimum value of \9\ at which the 
interface between the two phases may be quiescent). The halt of the front at 9 = —2 implies that the density behind 
it must take the conjugate value 9 = +2, which is indeed observed, see below. 

An intriguing case is \/3 < g < 2, when the nonsteady situation is possible whilst the maximum velocity of the 
steadily moving PT exceeds the sound velocity at 9 = 9*, suggesting that the PT front cannot send a sound wave 
ahead of itself. A typical picture for this case is displayed in Fig. 3, which shows that the acoustic precursor is 
nevertheless launched. The difference from the previous case is that the trough before the PT front quickly deepens, 
so that 9, rather than being stuck at the value quickly drops to the binodal value 9 — —2, see Fig. 2a. This 
leads to PT coming to halt earlier than in the case g < ^/3. Thus, an attempt to launch a "superfast" PT front into 
a deeply metastable phase produces an opposite result: the front generates strong acoustic waves and decelerates, 
seeing a less metastable state ahead of itself (physically, the metastability is relatively deep, as everything happens 
close to the critical point). 

An altogether different result is found at very small g, viz., onset of instability, which quickly switches the system 
into a fully turbulent state and destroys the interface. For instance, with the initial values 9q = —1-2 before the 
interface and 9 = 2.7 behind it, the instability sets in at g < 0.3. However, systematic study the present system as a 
model of ID turbulence is beyond the scope of this work. 

The situation behind the PT front deserves special description. We observe that ASW which hits the PT front in 
the last configuration shown in the upper part of Fig. 1 bounces from the front, travels to the right edge, is reflected 
from it, and then again hits the PT front (Fig. 4). A nontrivial feature of ASW in this case, evident in Fig. 4, is that it 
is reflected by the PT front as an antishock, or rarefaction wave. The present model does not admit steady rarefaction 
waves; however, the wave in Fig. 4 is unsteady, as it propagates between variable densities. Taking instantaneous 
values of the densities, we conclude that the antishock's velocity also obeys Eq. 

Analysis of the numerical results leads to the following inferences concerning the shuttling shock/antishock wave 
behind the PT front in Figs. 1 and 4: (i) bounces from both the right edge of the system and PT front do not change 
the absolute value of the velocity; (ii) bounces from the right edge are elastic without a change in the step height 
(9f — 9i); (iii) bounces from the PT front are highly inelastic: the shock bounces as an antishock and vice versa, each 
reflection decreasing the step height by a factor ~ 5, which explains why the shuttling wave disappears after a few 
reflections. 

Thus, we conclude that the system eventually drives itself into an equilibrium state with a quiescent interface 
between the liquid and vapor phases. A mechanism establishing the equilibrium is the generation of acoustic shocks, 
which bounce elastically from the system edge, but are muffled by strongly inelastic collisions with the interface. 
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FIGURE CAPTIONS 



Fig. 1. Propagation of the nonsteady condensation front with g = \/3. This and other figures display results for a 
typical case with the initial densities before and behind the front 9o = —1.2 and 82 = +2.7 (note that 9*{g = y/3) = 

— 1.63). Profiles 9(z) are shown through time intervals At, different in the upper and lower portions. An interval 

— 1 < 9 < 2.5, where nothing happens, is dropped. The bold arrows in the upper portion (and in Fig. 4) show 
velocities of the shuttling acoustic shock behind the front. 

Fig. 2. (a) The difference between the density 9, established ahead of the phase transition front after the passage 
of the acoustic precursor, and the critical density 9*. (b) The velocities of the acoustic shock and condensation front 
for g 2 = 2.4, 3.0, and 3.6. A small jump at t <~ 200 is due to hitting the front by the acoustic shock reflected from 
the right edge. 

Fig. 3. Evolution of the nonsteady condensation front at g 2 = 3.6. 

Fig. 4. The propagation of the antishock at g — y/3 between two inelastic collisions with the condensation front: 
before (a) and after (b) elastic reflection from the right edge. 
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